L  AD-R138  077 
UNCLASSIFIED 


INTEGRAL  GREEN'S  FUNCTION  SOLUTIONS  TO  DIFFERENTIAL 
EQUATIONS^)  AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON 
RFB  OH  SCHOOL  OF  ENGINEERING  U  H  AMELING  05  DEC  83 
AFI T/GEP/RH/83D-1  F/G  12/1 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL.  BUREAU  OF  STANDARDS-  19S3-A 


TIC  FILE  COPY 


AFIT/GEP/PH/83D-1 


INTEGRAL  GREEN’S  FUNCTION  SOLUTIONS 
TO  DIFFERENTIAL  EQUATIONS 

William  H.  Ameling 
Captain,  USAF 

AFIT/GEP/PH/83D-1 


Approved  for  public  release;  distribution  unlimited 


AFIT/GEP/PH/83D-1 


NUMERICAL  GRF.EN'S  FUNCTION  SOLUTIONS 
TO  DIFFERENTIAL  EQUATIONS 


THESIS 

Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science 

by 

William  H.  Ameling,  B.S. 

Captain,  USAF 

Graduate  Engineering  Physics 
December  1983 


■*?.  r  1  but  1  t  / 
Avnilut '  ;  1 .  i  v  •'  ■ 

A  v  >  \  '  •  .  r 


s 


g 


V  PREFACE 

/Ae  OLu-t^ors 

This  report  is  the  result  of  -ay  investigations  into  the 
numerical  solution  of  ordinary  and  partial  differential 
equations.  Green's  Functions  were  used  to  convert  the 
differential  equations  into  integral  form.  The  method  of 
central  finite  differences  is  a  common,  well  known 
method.^he  Method  of  Weighted  Residuals  was  used  as  a  third 
method  of  solution.  The  main  purpose  of  this  thesis  was  to 
compare  the  accuracy  and  advantages  and  disadvantages  of 
the  integral  equation  Green's  Function  method  to  the  other 

/At  iXsj.tl-*rs 

methods.  A  second  purpose  was  to  add  to  *ry  own  knowledge  on 

how  to  solve  differential  equations.  This  effort  also  served 

hn r  .  f/c, 

to  strengthen  -ay  own  abilities  to  write  computer  programs.  >1 

/■  '  -i-C  Aw  S 

-h**ePtried  to  give  enough  details  on  my  derivations  and 

equations  used  i*n  programing  that  others  could  use  them  in 

solving  their  own  problems  with  these  methods.  -  — - 

I  would  like  to  thank  Dr  Bernard  Kaplan,  my  adviser, 
for  his  guidance  and  advice  throughout  my  effort.  This 
thesis  was  sponsored  by  Nick  Pagano,  AFML/MLBM. 
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ABSTRACT 


Several  methods  of  numerically  solving  differential 
equations  in  one  and  two  dimensions  were  compared.  The  main 
method  used  was  an  integral  equation  solution  using  Green’s 
Functions  which  were  turned  into  matrix  equations  by  using 
the  Trapezoid  Rule.  The  method  of  finite  differences  was 
used  to  turn  differential  equations  into  matrix  equations 
which  were  solved  using  the  Gaussian  elimination  method. 

Both  a  Laplacian  Green’s  function  and  a  Helmholtz  Green's 
function  were  used  to  solve  the  one  dimensional  problem  by 
turning  it  into  integral  equations.  These  integral  equations 
were  turned  into  matrix  equations  using  the  Trapezoid  Rule. 
The  finite  difference  method  and  the  Laplacian  Green's 
function  method  gave  exactly  the  same  results.  The  Helmholtz 
Green's  Function  gave  slightly  better  results.  Two  types  of 
the  Method  of  Weighted  Residuals  were  studied:  the  Galerkin 
Method  and  the  Collocation  Method.  Both  of  them  gave  much 
better  results  than  finite  differences  and  the  Green's 
Function  methods  did.  The  Collocation  Method  gave  better 
results  than  the  Galerkin  method. 

For  the  two  dimensional  problem  an  attempt  was  made  to 
use  Method  of  Weighted  Residuals  to  reduce  the  partial 
differential  equation  to  an  ordinary  differential  equation 
which  was  then  salved  using  the  Green's  Function  method. 

Two  different  sries  were  used.  In  both  cases  the  results 

v  i 


were  very  bad  and  did  not  get  better  with  more  points.  The 
method  finite  differences  gave  reasonably  good  results  for 
the  partial  differential  equation.  The  method  of  separation 
of  variables  was  examined  to  see  if  it  could  be  used  to 
reduce  the  partial  differential  equation  to  an  ordinary 
differential  equation  whch  would  then  be  solved  by  the 
Green's  Function  Method.  The  analysis  showed  that  this 
approach  would  not  work  very  easily. 

The  programs  were  run  on  an  Apple  III  personal  computer 
and  were  written  in  UCSD  PASCAL. 
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1.  INTRODUCTION 


Many  of  the  real  world  problems  of  interest  to 
engineers  and  physicists  can  be  expressed  in  mathematical 
form  as  differential  equations  or  partial  differential 
equations.  These  equations  relate  differentials  of  functions 
in  time  and  space  to  each  other,  to  other  driving  forces  and 
to  boundary  or  initial  conditions.  This  thesis  is  only 
concerned  with  boundary  condition  type  problems. 

Boundary  conditions  can  be  classified  into  several 
types.  Dirichlet  boundary  conditions  specify  the  value  of 
the  unknown  function  everywhere  on  some  boundary.  Neumann 
boundary  conditions  specify  the  normal  derivative  (gradient) 
everywhere  on  the  boundary.  Mixed  boundary  conditions 
contain  Dirichlet  boundary  conditions  on  part  of  the 
boundary  and  Neumann  boundary  conditions  on  part  of  the 
boundary.  Cauchy  boundary  conditions  specify  both  the 
function  and  its  derivative  everywhere  on  the  boundary.  Only 
certain  types  of  boundary  conditions  will  lead  to  unique 
solutions  for  different  differential  equations  (Ref  5:706). 
Even  when  a  unique  solution  can  be  shown  to  exist,  finding 
it  may  be  very  difficult  analytically.  For  these  problems 
numerical  approximation  techniques  are  necessary. 

This  thesis  will  use  the  following  numerical  methods  to 


solve  differential  equations: 

(1)  The  method  of  central  finite  d if f erences(CFD) 

(2)  Green's  Functions 

(3)  The  Method  of  Weighted  Residuals 

Also  this  thesis  will  use  the  Gauss  Elimination  method 
to  solve  the  matrix  equations  that  arise  and  the  Trapezoid 
rule  to  numerically  evaluate  any  integrals. 

1.2  One-Dimensional  Case 

A  simple  ordinary  differential  equation  of  the 
Helmholtz  type  was  chosen.  This  permitted  the  use  of  two 
different  Green's  Functions  (Laplacian  and  Helmholtz),  the 
method  of  finite  differences  and  two  different  weighted 
residual  methods  (Galerkin  and  Collocation).  The  relative 
accuracy  of  each  of  these  methods  will  be  compared.  The 
Green's  Function  approach  relies  on  defining  a  special 
function  that  incorpates  the  boundary  conditions  of  the 
problem  and  which  permits  changing  the  differential  equation 
into  an  integral  equation  which  is  sometimes  easier  to 
solve.  The  method  of  differences  replaces  the  derivatives 
with  differences  between  the  values  of  the  unknown  function 
at  nearby  points  and  then  solves  the  resulting  set  of  linear 
equations.  The  Weighted  Residual  methods  approximate  the 
unknown  function  by  some  series  of  functions  with  unknown 
coefficients,  substitutes  it  in  the  equation,  and  applies  a 


weighting  factor,  and  then  tries  to  solve  the  resulting 
equations  for  the  coefficients. 

1.3  Two-Dimensional  Problem 


A  two  dimensional  heat  conduction  type  problem  was 
chosen  with  a  separable  driving  function  and  functional 
dependence  on  one  boundary  and  zeroes  on  the  other 
boundaries.  The  central  finite  difference  method  was  used 
again.  Then  attempts  were  made  to  reduce  the  partial 
differential  equation  to  ordinary  differential  equations 
that  could  be  solved  by  the  Green's  Function  Method. 

1.4  Purpose 


One  purpose  of  this  thesis  was  to  learn  more  about  some 
of  the  various  numerical  techniques  used  to  solve 
differential  equations,  do  numeical  integration,  and  solve 
large  sets  of  linear  equations.  Another  was  to  compare 
various  methods  to  each  other  in  hopes  of  learning  if  the 
Green's  Function  integral  methou  was  any  better  to  use  than 
the  CFD  method  which  is  in  common  use.  Hajdin  and 
Krajcinovic  (Ref  1  and  2)  contend  that  since  numerical 
integration  formulas  are  more  accurate  than  numerical 
d i f f er n t ia t ion  formulas,  that  methods  like  the  Green's 
Function  method  should  be  more  accurate  and  more  useful  in 
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complex  problems. 


1.5  Plan  of  Development 

First  the  theory  and  formulas  will  be  developed  for  the 
various  methods  for  the  one-dimensional  case.  Next  the 
algorythyms  and  computer  programs  will  be  written  and 
tested.  The  programs  will  generate  the  data  to  compare  the 
various  methods.  The  procedures  developed  here  will  be 
needed  for  the  two-dimensional  problem.  Then  the  theory  and 
formulas  for  the  two-dimensional  problem  will  be  developed. 
Finally  the  programs  for  the  two-dimensional  problem  will  be 
written,  tested,  and  used  to  generate  data. 

1.6  Equipment 

The  programs  will  be  written  in  PASCAL  and  run  on  my 
Apple  III  personel  computer  which  has  128K  bytes  of  memory. 


2. THEORY 


2.1  Green's  Functions 


Green's  Functions  are  a  common  tool  that  arise  in  many 
branches  of  physics.  Discussions  on  their  use  and  derivation 
can  be  found  in  most  math  physics  and  electromagnetism 
textbooks  (Ref  4-7).  They  are  found  in  the  solution  of 
differential  equations  and  integral  equations. 

One  way  of  deriving  them  is  by  making  use  of  the  Dirac 
Delta  function  &(<x).  The  Dirac  Delta  function  can  be 
represented  by  the  following  limit: 


86*0  - 


|  i  m 


Which  has  the  following  properties 


su)  =  o  *or<x.*o 

■*  00  ■f*/  ~  o 


(2.1.1) 


(2.1.2) 

(2.1.3) 


These  equations  lead  to  the  following  useful  result: 


-OQ 


(2.1.4) 


This  result  will  be  used  in  the  solution  of  the  following 


equation  by  the  method  of  Green's  Functions: 


asfr  =/"=  -f^y) 


(2.1.5) 


Boundary  Conditions:  y (o)  ci  y  yCO-b  (2.1.6) 


The  Green's  Function  G(x,x')  is  defined  by  the 
following  equations: 


of*,* 


-  Qf'1  -  -  $(**.-**') 


(2.1.7) 


Boundary  Conditions:  G(0.>x')  =  =  C>  (2.1.8) 


Multiply  Eq  2.1.5  by  G( x , x ’ )  and  Eq  2.1.7  by  y  (  x  )  and 
subtract  Eq  2.1.7  from  Eq  2.1.5  to  get 


6  y"-yG"  -  +  yMsi*-**)  (2.1.9) 


However 


&[*/“  y6'j=  s/-y&“ 


(2.1. 10) 


There  f  ore 


+  yMS(*-s. tO 


(2.1.11) 


*.:•%> 


Multiply  by  dx,  integrate  from  x=0  to  x=l,  and  rearrange 
terms  to  get 


yC*0*  |>y'- y&T  +  J  d* 


O  J  Q 


(2.1.12) 


Sustitute  the  boundary  conditions  (Eq  2.1.6  and  2.1.8)  to 


yC)t!)  *  cn  G 1  (  0^)  -  b  (a  ( f  o  ^  ) 
t  §  F(s#9y)  ct'Z 


(2.1.13) 


Where 


(2.1.14) 


A  more  general  result  for  three  dimensions  can  be  found 
by  making  use  of  Green's  Theorem: 

(  y  7*6  -  6  V*y)=  U  Cy  &  _  6£*J  (2.1.15) 


Where  n  is  the  normal  to  the  surface  S 


=  -  F(^sy) 

?*&  =  -  Z&-Z") 

for  x  on  the  surface  S  G  — 


(2.1.16) 

(2.1.17) 

(2.1. 18) 
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Use  Eq .  2.1.15  to  2.1.18  to  get  Eq  2.1.19  which  is  the  three 
dimensional  equivalent  of  Eq  2.1.12. 
y {$•)  =  -  (***£>  yU) 

+  J  ^  V  F(£,y) 


(2.1.19) 


2.2  Finding  Green's  Functions 


The  method  of  finding  the  Green’s  Function  G(x,x')  that 
satisfies  Eq  2.1.7  is  fairly  straightforward.  Let 


(.  ) 


•f  OV  >At/ 


(2.2.1) 


These  functions  satisfy 

the 

following 

equations: 

fGj  -o 

4 W 

(2.2.2) 

a**- 

(2.2.3) 

Which  have  the  general 

solut ions  of : 

6;  =  A  +  ior  <*><**/ 

G>»  =  C  +  D<* 


(2.2.4) 

(2.2.5) 


A,B,C,  and  D  can  be  found  by  applying  the  following 
conditions : 


,4 


w^m: 
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G,(0^«!)  =  I  >■*■')  =  o 

<B'/)  - 

^  6.(~,*0l  4- 6,  <*»•*'>  /  <  = 


-/ 


(2.2.6) 

(2.2.7) 

(2.2.8) 


Eq  2.2.7  is  true  because  we  want  the  function  G(x,x')  to  be 
continuous  over  the  entire  interval  x=0  to  x=l.  Eq  2.2.8 
comes  from  integrating  Eq  2.1.7  once  over  x.  Eq  2.2.8  shows 
that  the  derivative  of  G(x,x')  is  discontinuous  at  x=x'. 
Substitute  the  boundary  conditions  Eq  2.2.6  into  Eq  2.2.4 
and  Eq  2.2.5  to  get : 


f\  -  O  Ct  D  =  o 


(2.2.9) 


Use  Eq  2.2.7  in  Eq  2.2.4  and  Eq  2.2.5  to  get: 

B^e,;  =  C  t  OK  sfr'  ( D)  =  C  (2.2.10) 


Substitute  Eq  2.2.4  and  Eq  2.2.5  into  Eq  2.2.8  to  get: 


0-6  =  / 


(2.2.11) 


Eq  2.2.9,2.2.10,  and  2.2.11  have  the  solution: 


C  =  V  ,  D-  ,  6=  /-V 


(2.2.12) 


Which  yields  the  following  Green's  Function  solution  to  Eq 


2.1.7: 


L  0*  ( I  -4*)  «f 


C 1  -*J)  $  or 


a  /  ^ 


/  (2.2.13) 


Eq  2.2.13  will  be  referred  to  in  the  rest  of  thesis  as  the 
Laplacian  Green's  Function  since  it  is  the  solution  of  the 
one  dimensional  Laplacian  Equation  Eq  2.1.7. 

Another  Green’s  Function  that  will  be  useful  is  the 
solution  to  the  Helmholtz  Equation: 


Ct-4L> 


(2.2.14) 


With  the  same  homogenous  boundary  conditions  Eq  2.2.6.  This 
equation  has  the  general  solutions  of: 


b,  ~  A  sin^  t  8C05A  W  ^  <*J  (2.2.15) 

6*  s  c  t  Dc«5  4>  •f**'  (2.2.16) 


Substitute  Eq  2.2.6  in  Eq  2.2.15  and  Eq  2.2.16  to  get 


A  •  0  +  B  *  /  =  0  lJ-=0  (2.2.17) 

C5  i»  /  +Vco$l  =  0  D--Ctanl  (2.2.18) 

Substitute  Eq  2.2.15  and  Eq  2.2.16  into  Eq  2.2.7  and  Eq2.2.8 
to  get: 


A  =  C  £  /  -  £*+  fJ 
C  [c<>5*/+  S'*-*/  t«<*  f J-  ^  cold *-/ 


(2.2.19) 

(2.2.20) 


Eq  2.2.18,  2.2.19,  and  2.2.20  can  be  solved  to  give: 

«  /  .  y 

r;  5  i  1  , 

£05/*  TZlJ'"  (2.2.21) 

Eq  2.2.21  will  be  referred  to  in  the  rest  of  this  thesis  as 
the  Helmholtz  Green's  Function. 

A  useful  property  and  check  on  the  solution  is  that 
most  Green's  Functions  are  symmetric,  i.e.  G( x , x '  )  -G ( x '  ,  x  )  , 
when  the  differential  equations  and  boundary  conditions  are 
adjoint  (Ref  5:873-874). 


&  (*,*’)  =  [ 


2.3  Method  of  Weighted  Residuals 

A  second  method  of  solving  differential  equations  is 
the  Method  of  Weighted  Residuals.  Consider  the  following 
differential  equation  where  L  is  a  differential  operator: 

L  yM  t  F  (*)  =  0  (2.3.1) 

Expand  y  over  some  set  of  basis  functions  ^'f^to  get: 

N 

y  n  (d)  -  *2.  C;  (d) 


(2.3.2) 


fW9JW»vajr  3?  TF  TH  ■Tl’.  3  'A  LnW.rWTOrr,^.' 


Define  the  residual  R(x)  by: 

R  ('*)  £  L  t  F L*)  *  L  (  ^  C;  <£*  f  F£^t)  (2.3.3) 

*>*tci 

Let  V^)be  some^weighting  function.  Then  use  the  folowing 
set  of  equations  to  solve  for  Cj  : 

$'wj<*)K(*)d**0  *!,»>— ,#  <2-3-4> 

This  is  the  general  form  of  the  Method  of  Weighted 
Residuals.  There  are  various  names  given  to  different 
choices  of  weighting  functions  : 

Collocation  Method:  W-  (sft)  -  £  (2.3.5) 

j  *  y  J  J 

Galerkin  Method  Wj  (.*.)  *  rfj  ('&)  (2.3.6) 


2.4  Method  of  Finite  Differences 

A  third  method  to  solve  differential  equations  is  the 
Method  of  Finite  Differences.  This  method  approximates 
derivatives  with  differences  between  values  at  nearby  nodal 
points.  To  get  these  apr roximations  we  make  use  of  the 
Taylor  series  to  expand  y(x)  at  the  points  x-h,  x+h ,  and 

« 
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^LL 


x  +  2h : 


yu-h)  =yw->> y'U) * t y *U) ~^y"M (1 4  n 
+  -£?  y ^W-j£  Ao  +  oth*> 

y(-*+h)  =  yt*)  +  fey'i'-*)-*  £y'M  *  b?  y^c*)  (2.4.2) 

ti,  yCTM  +  Jj£  y*u>  +  0  Cb‘) 

yUtft,)  =  yf.0-*  zh'/'t*')  +  ai»*y%.)+ ^  V^y^u.")  (2.4.3) 

+  |  *>*yJP(*)+  ■&  ViY(»)  +  0( h‘) 

Where  0(h*)  refers  to  terms  containing  sixth  and  higher 
order  powers  of  h.  Now  subtract  Eq  2.4.1  from  Eq  2.4.2, 
subtract  Eq  2.4.2  from  Eq  2.4.3,  and  add  Eq  2.4.1  to  Eq 
2.4.2  to  get  the  following  equations: 

ihy'ud+^y**  W)  (2.4.4) 

lh)-yfc»>*H)=  k y'0>f  £ +0(h*)  (2.4.5) 

W*»+h)-h  =  Zyfa)  +  \^yllM  t  yi  ym(*)  (2,4*6) 

+  0(kfi) 

These  equations  can  be  solved  for  y’(x)  and  y''(x)  to  get: 


yU -  -£y"Mro(u')  (2 

y'4t)  =  r  <«♦«!- - 1  fc  y%)t  ef fc»j  <  ^ 


(2.4.7) 


.4.8) 


yV)  =  y  *)•»  jif  r  n  a  1) 

Vi*  lb 


(2.4.9) 


The  method  of  finite  differences  takes  h  to  bo  some  small 
number  and  neglects  the  second  terms  in  these  equations  to 
get  the  following  equations: 


1  3 


/fit-  h. 


(2.4.10) 

(2.4.11) 

(2.4.12) 


yfM~ 

y'M  = 
y*U)  & 


y(  _ 

v6**iQ 


* 


hx 


Eq  2.4.10  is  known  as  the  central  finite  difference  (CFD) 
approximation  for  y'(x).  It  has  an  associated  error  oi  oraei 

A, 

h  .  Eq  2.4.11  is  known  as  the  forward  finite  difference 
(FFD)  approximation  for  y'(x).  It  has  an  associated  error  of 
order  h.  Therefore,  using  the  CFD  method  with  symmetric 
points  around  x  should  give  more  accurate  results  than  the 
FFD  method.  Eq  2.4.12  is  the  CFD  approximation  to  y'*(x).  It 
has  an  associated  error  of  order  h  . 


2.5  Integration  Methods 


When  the  Green's  Function  Method  and  the  Method  of 
Weighted  Residuals  are  used  to  solved  differential 
equations,  integrals  will  need  to  be  evaluated  numerically. 
There  are  many  methods  of  numerically  calculating  integrals. 
Generally  they  work  by  taking  the  values  of  the  integrand 
f(x)  at  n+1  points,  finding  a  polynomial  of  nth  order  that 
fits  these  n+1  points  and  then  calculating  the  area  under 
this  polynomial.  This  procedure  is  repeated  for  small 
intervals  until  the  range  of  integration  is  covered. 


The  Trapezoid  Rule  is  the  simplest  method  of  doing 
numerical  integration  and  is  the  method  which  will  be  used 


in  this  thesis.  It  uses  first  order  polynomials  (i.e. 
straight  lines)  to  connect  points  on  the  curve  y=f(x).  These 
lines  form  trapezoids  whose  areas  are  easily  calculated.  The 
area  under  a  typical  trapezoid  with  a  base  length  of  h  and 
sides  of  f(xj,j  )  and  f(x^  )  is: 

Area*  (2.5.1, 


So  the  combined  total  area  of  all  trapezoids,  over  the 
interval  a  to^afeis: 

T«ia/  +  5  2) 

Therefore,  the  total  area  under  the  Trapezoid  Rule  is 


(2.5.3) 


yi  £[*(+,)+ &(*£*  —  +  *4  (<*n-i) 

*  ti  Ov)J 

This  approximation  has  an  associated  error  of  order  h3 . 
Simpson  *  s  Rule 

The  next  simplest  method  of  numerical  integration  is 
Simpson's  Rule  which  uses  polynomials  of  order  2  (i.e. 
parabolas)  to  fit  f(x).  It  can  omly  be  used  the  number  of 
intervals  used  to  divide  x=a  to  x=b  is  even.  The  resulting 
equation  is  (Ref  6:350-351): 

jbfuv*  *  \  [-f c*. «(•,)+*«*,) 

a  *  •••  +  <2-5-4) 

This  approximation  has  an  associated  error  of  order  h^  .  A 
earlier  thesis  (Ref  3:26-30)  found  that  the  trapezoid  rule 
gave  more  accurate  results  than  Simpson's  Rule  for  problems 
involving  linear  Green's  functions.  This  is  because  at  the 
point  x=x'  a  parabola  includes  more  incorrect  area  than 
trapezoids  do  for  a  linear  Green's  Function.  Therefore, 
because  it  is  simpler,  can  be  used  for  any  number  of  nodal 
points  (Simpson's  Rule  requires  even  number  of  points),  and 
because  it  gives  better  results  for  the  Laplacian  Green's 
Function  this  thesis  will  use  the  Trapezoid  Rule  when  doing 
numerical  integration. 


In  sections  2. 1,2. 4,  and  2.5  we  have  discussed  various 


methods  of  solving  differential  equations.  It  will  simplify 
notation  and  make  it  easier  to  see  the  methods  of  solution 
by  adopting  matrix  notation  and  methods.  Consider  the 
following  equation: 

y  (&)-  F(«)t  §  Cjfc)*/*/  (2.6.1) 


Use  the  Trapezoid  Rule  Eq  2.5.3  to  replace  the  definite 
integral : 

y(*i)*-F6w)+  i  Dj (2  6  2) 

D,=  =  i  0*»  »  Dv-I  *  k  (2.6.3) 

Irt  ylf.i')  i  y; ;  Ff^o*  (j  j  s  G;j 

■f  (-*>•)  a  t-' 

Use  these  results  and  definitions  to  write  the  following 
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matrix  equation: 


&// 4/ 


&vrfj 


•  G/A^A/ 


(2.6.5) 


In  more  compact  form  Eq  2.6.6  can  be  written  as: 


y  =  F  +■  D 


(2.6.6) 


Another  equation  that  will  be  needed  is  the  following 
integral  equation: 


y  (*s)  =•  F(^k)  t  j 


(2.6.7) 


Using  the  Trapezoid  Rule  Eq  2.6.7  can  be  written  as: 


y<  =  F;  +  £  G*  Oj  yj 


(2.6.8) 


Or  in  matrix  form  as: 


y  =  F  +  G  J)  y 


(2.6.9) 


Where :  ^  — 


Sa/|  -  •  •  G^A/ 


(2.6. 10) 


p  -  ( V  0  \B  0,  s(J7 

V  0  'oN 


(2.6.11) 


Eq  2.6.9  can  be  rewritten  as: 


[l-G  0]y  =  F 


(2.6.12) 


2.7  Converting  CFD  Into  Matrix  Eouat i( 


In  section  2.4  the  formulas  for  central  finite 
differences  were  developed.  In  this  section  they  will  be 
used  to  turn  the  following  differential  equation  into  a 
matrix  equation: 


+  g&O  yfa)  =  f  (#) 


(2.7.1) 


Sustitute  Eq  2.4.12  into  Eq  2.7.1  to  get: 


,  3G) yU)  -  «*) 


(2.7.2) 


Make  the  following  notation  changes: 
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yk*)?  /•  5  5  f  ^')-£  -f; 

(2.7.3) 


The  resulting  general  equation  is: 

y,>i  +y;-i  +  b  ■f*  c 2.7.4) 

■far  1  »  3,  '  /V"/ 

Let  the  boundary  conditions  on  y(x)  be: 


y  1  ^  A  yv  5  E> 

(2.7.5) 

Eq  2.7.4  and  Eq  2.7.5  generate  the  following  set 

of  linear 

equations : 

y. 

~  A 

yf  t 

= 

y»  r  Yi 

= 

• 

«  # 

• 

« 

(2.7.6) 

1 

«* 

_c 

II 

£ 

•*» 

yv  «  P 


Eq  2.7.6  can  be  written  in  either  of  the  following  matrix 
forms : 


10  0  O  o 

I  /  °  0  *  *  ‘ 

0  I  ($yt)  10  •  •  • 

* 

•  i  1 

,  O  O  I 


A  y  =  b  (2.7.9) 


2.8  Matrix  Inversion 

In  sections  2.6  and  2.7  the  following  form  of  matrix 
equation  has  arisen  (Eq  2.6.12  and  Eq  2.7.9): 

A  *  -  &  (2.8.1) 

Where  A  is  the  coefficient  matrix,  x  is  the  unknown  column 
matrix,  and  B  is  constant  column  matrix.  The  problem  is  to 
determine  x. 

The  most  commonly  employed  method  is  the  Gauss 


elimination  method  (Ref  4:1-4).  This  method  makes  use  of  the 


fact  that  the  following  algeraic  manipulations  performed  on 
a  set  of  linear  equations  leaves  their  solution  unchanged: 

(1)  Multiplication  or  division  of  any  equation 
(row)  by  a  constant 

(2)  Adding  or  subtracting  one  equat ion ( ro w )  to 
another  equation(row). 

Consider  as  an  example  the  following  set  of  four  equations: 


An  4}3l  Aft 

Am  A**  Am 

Ay  A-J3  A** 
Am  Ah 


(2.8.2) 


The  Gauss  elimination  method  is  simply  a  sequential 
application  of  row  operations  (1)  and  (2)  above.  The  goal  is 
to  reduce  all  elements  below  the  diagonal  elements  to  zero. 
First  subtract  A>i/A,i  times  the  first  row  from  the  second 
row,  subtract  Av/A//  times  the  first  row  from  the  third  row, 
and  subtract  Ami/A  „  times  the  first  row  from  the  fourth  row. 
This  produces  Eq  2.8.3  where  all  the  elements  in  the  first 
column  are  zero  except  for  the  first  element: 

An  A Aft  A/*j\  /^A 

0  An  A  *3  A^tf 

0  AVa  A '„A'„  5:1  ft,'  (2.8.3) 

V  A 4 A  Any  Ahh  W*  \  fy/ 


V.V.  . 


/• 


Next  start  with  the  second  row  and  follow  the  same  steps  to 
reduce  the  third  and  fourth  elements  of  the  second  column  to 
zero.  Finally,  use  the  third  row  to  reduce  the  fourth 
element  of  the  third  column  to  zero: 


%|  A|3 

0  0  Aji 

\o  o  o 


(2.8.4) 


Eq  2.8.4  can  now  be  used  to  find  the  solution  matrix  by 
simple  substitution: 

^  b"/C ,  ,  (2.8.5) 

B3  -3  (2.8.6) 

(  6#.  ~  A*/*5*  -  (2.8.7) 

(2  8  8  ) 

"^7  =  C  B/  ■  Aii  s#*  -  ~  An i 

During  these  operations  the  key  elements  are  the  diagonal 
elements  which  are  known  as  the  pivot  elements.  The  method 
does  not  work  if  a  pivot  element  is/becomes  zero.  Also  if 
the  pivot  element  is  small  compared  to  the  other  elements 
below  it  it  will  led  to  round-off  errors.  The  solution  is  to 
switch  rows  around  to  put  the  largest  element  in  the  column 
at  or  below  the  diagonal  in  the  diagonal  position.  This 


comparison  is  done  when  the  method  is  done  with  the  the 
previous  column. 

Another  thing  to  note  is  the  types  of  matrices  to  be 
solved  are  different.  The  Green's  Function  meLhod  and  the 
Method  of  Weighted  Residuals  generate  'dense'  matrices  where 
most  of  the  elements  are  nonzero.  The  central  finite 
difference  method  on  the  other  hand  generates  a  tridiagonal 
matrix  where  only  the  diagonal  elements  and  the  elements  to 
each  side  are  nonzero.  The  solution  of  the  tridiagonal 
matrix  is  much  simpler  because  of  the  fewer  number  of 
operations  needed  (of  order  N  versus  order  N^)  and  computer 
memory  storage  requirements  can  be  reduced  drastically  from 
N*N  to  3N.  Also  because  of  the  fewer  number  of  operations 
needed  to  solve  a  tridiagonal  matrix,  the  accumlated 
round-off  errors  will  be  much  smaller. 


3.  ONE  DIMENSIONAL  PROBLEM 


3.1  Computer  Programs 


The  computer  programs  used  to  solve  a  one  dimensional 
problem  for  the  various  methods  used  were  all  very  similiar 
The  flow  diagrams  a.ppears  in  Figure  3. 


Read  N:  number  of  data  points 
Read  Y0,Y1:  boundary  conditions 

Generate  array  of  nodal  points 
Generate  Green's  array 

Generate  numerical  integ.  coef  .  array 

Generate  the  coef.  array  and  the  constant 
array  using  matrix  operations 

Use  Gauss  Elimination  method 
procedure  to  solve  the  matrix  eq. 

Generate  exact  solution  array 
and  compare  to  get  error  array 

Print  out  solution  array  and  error  array 


Fig.  3  Flow  Diagram  for  Computer  Programs 


All  that  needed  to  be  changed  for  each  method  used  were 
the  procedures  that  generated  the  matrices  needed  and  the 
appropiate  products  of  matrices. 

3.2  Problem 


The  following  one  dimensional  differential  equation 
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will  be  solved  by  the  Green’s  Function  method  using  both 
Laplacian  and  Helmholtz  Green's  Functions,  the  method  of 
central  finite  differences  (CFD),  and  by  the  Method  of 
Weighted  Residuals  (MWR)  using  collocation  and  Galerkin: 

"U)=-yM+-*  <3.2.1) 


Boundary  Conditions: 


y(o)  -  yd)  -  O 


(3.2.2) 


This  equation  was  chosen  beause  it  is  the  simplest  second 
order  differential  equation  whose  Laplacian  Green's 
Function  solution  gives  a  integral  equation  and  also  has  a 
simple  Helmholtz  Green's  Function  solution.  It  will  also 
show  tfie  essential  features  of  CFD  and  MWR  methnods. 

The  solution  to  Eq  3.2.1  is  relatively  easy.  It  is  to 
see  that  the  particular  solution  y p(x)  =  x  satisfies  the 
inhomogeneous  Eq  3.2.1  This  leaves  the  solution  of  the 
following  homogenous  equation  to  be  found: 

y"  +  y  s  O  (3.2.3) 


Eq  3.2.3  can  be  solved  by  assuming  that  y(x)  is  a  linear 
combination  of  exponentials.  Substitute  the  following 
equation  (Eq  3.2.4)  in  Eq  3.2.3  to  get  Eq  3.2.3: 

V'X 


y  M  -  c 


(3.2.4) 


§ 


(D  *tl)eD~»0 


(3.2.5) 


Divide  by  the  exponential  to  get: 

0*t<  -O  «>  D  = 


(3.2.6) 


Therefore  the  homogenous  solution  is: 
y  H  i'X)  -  he.  frco^jt  f  (3.2.7) 


Therefore  the  general  solution  is: 


yU)-  y?t  Yh  - 


sx  +  A  -r  Bji’w/ar 


(3.2.8) 


Apply  the  boundary  conditions  Eq  3.2.2  to  get: 

y(0)~  0  -  0  +  A* /  + 
y  (i)  -  0  s  /  +  A  co*/  t  B  si*»  / 


Which  have  the  solution: 


A=  0  a*  c/  6  = 


5<n  / 


(3.2.11) 


Which  gives  the  general  solution  as: 


=  frrf 


(3.2.12) 


3.3  Green’s  Function  Integral  Equations 

We  know  from  section  2.1  and  Eq  2.1.5  and  Eq  2.1.12 
that  the  Eq  3.2.1  can  be  put  in  the  following  form: 

-  [Gy7-  J0  (3.3.D 

^  Cy  6  6*^0  <5^ 

Use  the  boundary  conditions: 

y&0~  y(<)  =  0 

'  (3.3.2) 


to  simplify  Eq  3.3.1  to  the  following  form: 

y  6#/)  -  ( y  -  G 


The  Laplacian  Green's  Function  is  (from  Eq  2.2.13): 


'ifc  ( I  'sat,1) 

<*' ( I  '<0 


•fck' 

fey 


(3.3.4) 


Eq  3.3.3  can  be  turned  into  the  following  matrix  equation 
(see  section  2.6): 


Eq  3.3.5  can  be  put  in  the  following  form  using  matrix 
methods : 


(x-  6|-D)y  •  G 


Eq  3.3.7  has  the  form  of: 

A  7  =  "5 


(3.3.8) 


Where 


/“»  -  J 


(3.3.9) 


Eq  3.3.8  can  be  solved  using  Gaussian  elimination 
(see  section  2.8)  to  get  the  Laplacian  Green's  Function 
solution  to  y(x). 

Another  method  of  solving  Eq  3.2.1  is  to  use  the 
Helmholtz  Green's  Function.  The  Helmholtz  Equation  is: 


(3.3.10) 


The  Helmholtz  Green's  Function  is  (Eq  2.2.21): 


<.-.7  .7  ..■  .7  .7  .7 -. 


r  5<h  yja  Cc  0$  yX>  **  >2-^—  J  ^OV  «*<\X' 

G*C*s+)-  7  *«  ' 

v  (  •  ,r  (3.3.11) 


The  equivalent  of  Eq  3.3.3  is  the  following  equation: 


rJ  ,  / 

y  ( <jb)  —  ~  J  'Z, 


(3.3.12) 


Which  has  the  matrix  form  of  (using  the  Trapezoid  Rule) 


=  -  &  D 


<3.3. 13) 


Eq  3.3.13  can  be  solved  easily  by  matrix  multiplications 
without  having  to  use  Gaussian  elimination. 


3.4  Central  Finite  Difference  Solution 


The  next  method  used  to  solve  Eq  3.2.1  is  the  method 
of  central  finite  differences  (CFD)  (see  sections  2.4  and 
2.7).  Use  Eq  2.4.12  to  get: 


y*=  y^JYHZhi  =  -y;  +  X; 


(3.4.1) 


This  equation  simplifies  to  the  following  set  of  linear 
equations: 


y,-_l  +  (k’-*)y;  + 


(3.4.2) 


With  the  boundary  condition  equations: 


y/  =  yv  -  ° 


(3.4.3) 


Eq  3.4.2  and  Eq  3.4.3  can  be  turned  into  the  following 
matrix  equation  and  definitions: 


(3.4.4) 


//  0  o  o  0 

I  l  I  o  0 

o  I  (w^> 1  0 


\ 


i  a.'-# 

CJ  o 


y  =  Cy< y*J 

B  -  [  X,  ■  •  ■  XtfJ 


(3.4.5) 


(3.4.6) 

(3.4.7) 


Eq  3.4.4  can  be  solved  by  the  Gauss  elimination  method  to 
get  the  solution  for  y(x)  at  the  nodal  points. 


3.5  Method  of  Weighted  Residuals  .Solution 


Eq  3.2.1  can  also  be  solved  using  the  Method  of 


.  \ 

v' 


V 


Weighted  Residuals  (see  section  2.3).  Chose  a  set  of  basis 
functions  which  satisfy  the  boundary  conditions 
y(0)=y(l) =0 : 


(3.5.1) 


Substitute  Eq  3.5.1  into  Eq  2.3.2  to  get  an  approximation 
to  y ( x ) : 


ym  M  *  Z  1 


(3.5.2) 


The  residual  is: 


f =  /»*.  +  Y**>  ** 


(3.5.3) 


t  ajCj(j*0*J'*'jQ^yi] 

*  t  %  (/-<&)  -  SZ 


(  3 . 5  .  A  ) 


y-i  + 


(3.5.5) 


Ga lerkin  Method 


For  the  Galerkin  Method  we  use  the  functions  in  Eq 
3.5.1  as  the  weighting  functions  and  integrate  from  x=0  to 
x  =  1  and  set  the  Residual  equal  to  zero  to  get: 
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-  -»  -*  J 


f  = 0 

J  n 


(3.5.6) 


Substitute  Eq  3.5.5  to  get 


i  r  ** 

o  =  W  *(!-■»)$-*’ + 1 

)o  -+0*01] 


(3.5.7) 


Simplify  Eq  3.5.7  to  get 


Use  the  following  result  (Eq  3.5.9)  to  simplify  Eq  3.5.8 


and  get  Eq  3.5.10: 


*>■«"  •  f *nl  =  -/ 


(3.5.9) 


|  a- -  3i!  +  (3.5. 1 0) 

~  TTfT^  *  i-*6i2>  3 

Eq  3.5.10  can  be  turned  into  the  following  matrix 
equation  and  definitions: 


fit  a-  B 


...  -  jCj-0.  t  _i 

A,j  "  iij-i  ;-»j  T  if j+  (  iijta 


6,  = 


v  2  «  ▲  ^ 


(3.5.11) 


(3.5.13) 

(3.5.13) 


Where  #  is  the  column  matrix  of  Ct^  .  Eq  3.5.11  can  be 
solved  using  the  Gauss  elimination  method. 

Collocation  Method 


The  collocation  method  uses  a  set  of  displaced  Dirac 
Delta  functions  as  the  weighting  functions  in  Eq  3.5.6: 

f 1 R*  M  -  0  1=1)2,--'  jA/ 

J O  (3.5.14) 


Substitute  Eq  3.5.5  into  Eq  3.5.14  to  get: 

0  =  +  fVW,)*  j  .(-,/)}]<  3. 5. 15) 

Eq  3.5.15  can  be  put  into  the  form  of  Eq  3.5.11  with  the 
following  definitions: 

A;j  s  X\*  ( I  -<*; )  +  j  J  fr-  /*'*>;  fjt  0 J 

6;  =•  t 

Eq  3.5.11  can  again  be  solved  by  the  Gauss  Elimination 
method  to  get  coef  Qj  . 

To  get  y(x)  for  both  Galerkin  and  Collocation  methods 


(3.5.16) 

(3.5. 17) 


substitute  the  solutions  into  Eq  3.5.2  and  the  values  of 
the  nodal  points  in  x: 


mssm  mssm 


v  v  V  v  .*  v  -*  .*  - 


m 

Y«l  (**•)“  Z 


3.6  Numerical  Results  and  Comi 


(3.5.18) 


This  section  will  present  the  numerical  results  and 

errors  of  the  Laplacian  Green's  Function,  Helmholtz 

Green's  Function,  central  finite  difference  (CFD), 

Galerkin,  and  Collocation  methods. 

First,  Table  1  shows  results  from  the  CFD  method,  the 

Laplacian  and  Helmholtz  Green's  Function  methods  and  the 

exact  value  of  y(x).  Results  are  shown  for  one,  three  and 

five  interior  nodal  points.  Second,  Table  2  shows  the 

relative  percentage  errors  of  the  results  in  Table  i. 

Table  1.  CFD,  Laplacian,  and  Helmholtz  Results  for 
Different  N  and  at  Several  Points 


Y  CFD 


Y  LAP 


Y  HELM 


Y  EXACT 


1  0.5000  -0.07412857  -0.07142857  -0.06828782  -0.06974697 


3  0.2500 
0.5000 
0.7500 


-0.04427401 

-0.07015590 

-0.06040304 


-0.04427401 

-0.07015590 

-0.06040305 


-0.04378418 

-0.06938332 

-0.05974305 


-0.04401368 

-0.06974697 

-0.06005621 


5  0.1667 
0.3333 
0.5000 
0.6667 
0.8333 


-0.03056336 

-0.05564811 

-0.06992782 

-0.06837619 

-0.04640671 


-0.03056336 

-0.05564810 

-0,06992781 

-0.06837618 

-0.04640670 


-0.03041289 

-0.05537471 

-0.06958544 

-0.06804304 

-0.04618201 


-0.03048348 

-0.05550328 

-0.06974697 

-0.06820101 

-0.04628921 
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Table  2.  Relative  Percentage  Errors  of  CFD, 
Laplacian,  and  Helmholtz  Methods  for 
Different  N  and  at  Several  Points 


X 

Y  CFD 

Y  LAP 

Y  HELM 

0 . 5000 

2.411 

2.411 

-2.092 

0.2500 

0.591 

0.591 

-0.521 

0.5000 

0.586 

0.586 

-0.521 

0.7500 

0.577 

0.577 

-0.521 

0.1667 

0.262 

0.262 

-0.232 

0.3333 

0.261 

0.261 

-0.231 

0.5000 

0.259 

0.259 

-0.232 

0.6667 

0.257 

0.257 

-0.232 

0.8333 

0.254 

0.254 

-0.232 

Table  3.  Relative  Percentage  Error  Trends  at  X=.5 
for  Different  N  for  CFD,  Laplacian,  and 
Helmholtz  Methods 


N 

H 

Y  CFD 

Y  LAP 

Y  HELM 

1 

0.50000 

2.41101 

2.41101 

-2.09206 

3 

0.25000 

0.58629 

0.58630 

-0.52138 

5 

0.16666 

0.25929 

0.25928 

-0.23159 

9 

0.10000 

0.09313 

0.93086 

-0.08335 

19 

0.05000 

0.02343 

0.02313 

-0.02094 

29 

0.03333 

0.01052 

0.01039 

-0.00930 

Table  3  shows  the  relative  errors  at  x=0.5  as  the  number 
of  nodal  points  (N)  increases  and  the  step  size  (H) 
decreases.  It  can  be  seen  that  the  size  of  the  errors 
decreases  for  smaller  steps  as  would  be  expected.  The  most 
unexpected  result  in  Tables  1,2,  and  3  is  that  the  CFD 
method  and  the  Laplacian  Green's  Function  method  give 
exactly  the  same  results  (differences  in  the  last  decimal 
place  are  due  to  different  accumulated  round-off  error).  A 


similiar  result  was  gotten  in  a  previous  thesis  (Ref 


3:62-65).  The  Helmholtz  Green's  Function  method  gave 
slightly  better  results  in  all  cases  than  the  CFD  and 
Laplacian  Green's  Function. 


Table  4.  Laplacian,  Galerkin,  and  Collocation  Results 
for  Different  N  and  at  Several  Points 


Y  LAP 


Y  GAL 


Y  COLLOC 


Y  EXACT 


3  0.2500  -0.04427401 
0.5000  -0.07015590 
0.7500  -0.06040304 

5  0.1667  -0.03056336 
0.3333  -0.05564811 
0.5000  -0.06992782 
0.6667  -0.06820101 
0.8333  -0.04640671 


-0.04403237 

-0.06974640 

-0.06003848 

-0.03048322 

-0.05550319 

-0.06974734 

-0.06820088 

-0.04628902 


-0.04396500 

-0.06974170 

-0.06009403 

-0 . 03048368 
-0.05550334 
-0.06974698 
-0.06820092 
-0.04628904 


-0.04401368 

-0.06974697 

-0.06005621 

-0.03048348 

-0.05550328 

-0.06974697 

-0.06820101 

-0.04628921 


Table  5.  Relative  Percentage  Errors  of  Lapacian, 
Galerkin,  and  Collocation  Methods  for 
Different  N  and  at  Several  Points 


m 


X 

Y  LAP 

Y  GAL 

Y  COL 

0.2500 

0.59147 

0.04247 

-0.11058 

0.5000 

0.58630 

-0.00081 

-0.00756 

0.7500 

0.57751 

-0.02952 

0.62973 

0.1667 

0.26204 

-0.00087 

0.00025 

0.3333 

0.26095 

-0.00015 

0.00011 

0.5000 

0.25929 

0.00053 

0.00001 

0.6667 

0.25687 

-0.00019 

-0.00013 

0.8333 

0.25385 

-0.00041 

-0.00036 
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Table  6.  Relative  Percentage  Error  Trends  at 
x=0.5  for  Different  N  for  Laplacian, 
Galerkin,  and  Collocation 


N 

H 

Y  LAP 

Y  GAL 

Y  C0LL0C 

3 

0.25000 

0. 58629 

-0.00081 

-0.00756 

5 

0.16666 

0.25929 

0.00053 

0.00001 

9 

0.10000 

0.09313 

0.00037 

-0.00001 

19 

0.05000 

0.02342 

-0.00064 

-0.00003 

29 

0.03333 

0.01052 

0.00004 

0.00002 

From  Table  5  it  can  be  seen  that  the  Galerkin  and 
Collocation  methods  are  much  more  accurate  then  the 
Laplacian  Green's  Function  Method  (and  CFD  method). 
However,  the  size  of  the  error  varies  by  a  large  amount  at 
various  nodal  points.  From  Table  6  it  can  be  seen  that 
there  is  nothing  to  be  gained  by  going  to  large  numbers  of 
nodal  points  since  the  size  of  the  error  stopped 
decreasing  at  N=9  for  the  Galerkin  method  and  at  N=5  for 
the  Collocation  method.  One  potential  problem  with  using 
the  Galerkin  method  in  more  complicated  probims  is  that 
doing  the  integral  in  Eq  3.5.6  will  be  more  difficult.  In 
fact,  numerical  integration  maybe  necessary  to  evaluate  Cq 
3.5.6  and  this  will  decrease  the  accuracy  of  the  final 
results.  On  the  other  hand,  the  Collocation  method  does 
not  suffer  from  this  problem  of  evaluating  Eq  3.5.6  since 
it  uses  Eq  3.5.14  instead  which  is  always  easy  to 
integrate  because  of  the  Dirac  Delta  function. 


4.  TWO  DIMENSIONAL  DIFFERENTIAL  EQUATION 


4 . 1  Problem 


The  following  form  of  equation  often  arises  in  heat 
conduction  problems  and  in  electromagnetism: 


V*T*K»T  =  ^+5;tKaT  »  * 


(4.1.1) 


If  k*k  takes  the  following  form  the  above  equation  can  be 
separated : 


+  $00 


(4.1.2) 


Let  T ('&,'/)  ~  *(+)y(y) 


(4.1.3) 


Susitute  Eq  4.1.2  and  Eq  4.1.3  into  Eq  4.1.1  to  get 


+  O*-  b*  +  3^)]  y fy)  *  o 


(4.1.4) 

(4.1.5) 


To  get  a  problem  that  would  result  in  an  integral  equation 
let  us  start  with  a  T(x,y)  and  find  out  what  k*k  is: 

n-f  (4.,.6) 

£1  =  *  cat*  (4.,.7) 


gy-  =  L*  W  T  “  J  J/  S'**.  (A. 1.8) 

Therefore: 

?!r,?!r  _  -*T  -  vy  *'•»*£  *’mtl 

'  *■  ga>  (4.1.9) 

7  tirycoi*? 

Use  Eq  4.1.10  and  Eq  4.1.11  to  simplify  Eq  4.1.9  into  Eq 


1.12: 

Try 

5'« 

=  U  r.« 

*L 

(4.1. 

.  10) 

7T**  C* 5 
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005 
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Eq  4.1.12  is  the  two  dimensioal  differential  equation  which 
will  be  used  to  try  out  the  methods  of  solving  it.  The 
boundary  conditions  are: 

T(0,y)=  TCI>y)=  =  O  (4.K13) 

()=  ^  (4.1.14) 

The  exact  solution  to  Eq  4.1.12  is  Eq  4.1.6.  Eq  4.1.12  can 
be  put  inthe  following  more  compact  form: 


F(*,y>TU,y)  =  o 

v  ^  T*  7T  ^  W  ZT^+ 


M  / 


(4.1.15) 

(4.1.16) 


+  Tit//  +  Tr/-/  +  Tr/t/^Tv 


(4.2.5) 


Multiply  by  h*h  to  get: 


Tr-i,/  ■*■  trfVfj j-*)Tv»o  (4.2.6) 


Eq  4.2.6  upon  careful  examination  can  be  converted  into  the 
following  matrix  equation: 


l  O  l  0  O  O  o  o 


I  l  O  l  0  o  O  o 

o 


\ 


0  I  (h%-*)  0  o  I  0  0 

o  O  O  I  o  0 

I  0  I  (h*F5r«)  I  0  I  0 
oioi  (PMo  o  i 
0  0  I  O  0  »-H)  I  0 

o  o  G  I  G  I  (k*F ^ 
0  0  0  0  J  0  l 


I 

6 

0 

c> 

0 

0 


M 

-To 

\  T*y 

l 

1 

T3a 

-t„ 

0 

TM 

-7 *rV 

v  I7''3  J 

-tJ3 

>/W 

~  T<<f"T5Y  j 
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The  ones  are  located  next  to  the  diagonal  and  at  N  elements 
away  from  the  diagonal  element.  Any  of  these  elements  will 
be  zero  instead  if  they  correspond  with  elements  on  the 
boundary  of  the  problem.  The  following  conditions  incorpate 
the  boundary  conditions  (see  Eq4.1.13  and  Eq  4.1.14): 


T) j  =  t5-«j  =•  >  -  O 

Tr*  =  *,c**(5£0 


(4.2.8) 

(4.2.9) 
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Eq  4.2.7  has  the  form  of: 

%  T  =  S' 


(4.2.10) 


Eq  4.2.10  can  be  solved  by  the  Gauss  elimination  to  get  the 
solution  for  T  at  the  nodal  points. 


4.3  Method  of  Weighted  Residuals  Solution 


Let  T(x,y)  be  approximated  by  the  following  expansion: 


1  u-i 


(4.3.1) 


Then  use  Eq  4.3.1  to  reduce  Eq  4.1.13  to  an  ordinary 
differential  equation  which  will  be  solved  by  the  method 
Green's  Functions: 


2  -*■'{/-■■*)  "(y) 

=  £  c.(y)  -.*!'*[  ;-i 


rJ 

t 

«/=/ 

The  residual  R(x)  is 


(4.3.2) 

(4.3.3) 


Ux,y-)  =  £&£  t  t  F(*>Y)  T^y) 

1  fiy*  / 


(4.3.4) 


Substitute  Eq  4.3.2  and  Eq.  4.3.3  into  Eq.4.3.4  and  simplify 


43 


to  get: 


RU>y)=  f.  +  C!(y)C' 

i9,L  *)]} 


3.5) 


Multiply  by  the  shifted  Dirac  Delta  weighting  functions  and 
integrate  to  get  (see  Eq  3.5.1  ): 


0  =  2  (4.3.6) 

+  c;(  y)[;^j  (!*0) 


The  only  way  that  this  equation  can  be  solved  is  if  each 
term  for  a  different  i  value  is  equal  to  0: 


+  •' ( »>0j rF(*i>y)'*5 


(4.3.7) 


This  equation  can  be  simplified  into  the  following 
differential  equation: 


C;"(f)  r  C;(y)  f,j(y)  -  C 


(4.3  O 


Where  f  ij  (y)  “  1  C  *  f 


(4.3.9) 


Next  solve  Eq  4.3.8  using  a  Laplacian  Green's  Function  to 
get  ( see  Eq  2.1.12): 


C\(y)=  C Gc/(y')-  C;(y)^r'&ly^a 
f\0  CfCy')  fa(/)  &Cy>y/)*ty/ 


(4.3.  10) 


Boundary  Conditions  are: 


&Cy>0  =  o  ■,  C;(o) » o 

C:  O)  y  O 


(4.3.11) 

(4.3.12) 


The  Laplacian  Green's  Function  is  (see  Eq  2.2.13): 


s,w)=  f  y,c';y,J  w 

L  Y  Cl- y)  W  y  >  y' 


y>r 


(4.3.13) 


The  derivative  of  the  Green's  Function  at  y'=l  is: 


dGJy,y')l  ^  -y 


(4.3.14) 


Therefore  Eq  4.3.10  becomes: 

C;Cy)  -  yC;(i)t  f'  C; (y)  -f  ;j (y) 6 (yt/)4{ 4.3.15, 


This  can  be  put  into  the  matrix  form: 

Zi  -  y  c,(t)  +■  (6f)  p  (4.3.16) 


Therefore: 


( J  -  (6f3  D  )  C;  =*  y  c-;(0 


(4.3.17) 


To  find  the  boundary  conditions  at  y=l  use  the  following 


5 

fr,  (4.3.18) 

This  equation  can  be  turned  into  the  following  matrix 
equation : 

A  c~,u)  ~  3 

;  -jr  (4.3.  19) 

A,0  =  Xj)  ^  J  (A.  3. 20) 

Eq  4.3.19  can  be  solved  by  Gauss  elimination  to  get  the 
to  substitute  into  Eq  4.3.17  which  can  then  be  solved  for 
.  The  main  probem  with  this  method  is  that  turns  out  to 
be  a  function  of  the  x  used  in  the  displaced  Dirac  Delta 
Functions.  As  a  result  we  have  to  solve  N  ordinary 
differential  equations  for  for  each  x  for  a  total  of  N%N 
differentials  equations  solved  (where  N  is  the  number  of 
nodes  along  each  axis). 

4.4  Alternate  MW  R  Solution 

Another  solution  to  Eq  4.1.15  uses  the  following 


expansion: 


aj  * 


expansion  for  T(x,y): 


(4.4.1) 


7 ~,j  (* ,y)  =  f  y*  ( ;  *•*)  c;  <V) 

i~  I 

This  expansion  gives  the  following  expansion  for  the 
residual  : 

fUx^y)  =  £  [c,"(y)  -  F6r,y)C, •(■/)]  (4.4.2) 

,>l  •  5**  (.iff*#') 

Use  the  displaced  Dirac  Delta  fuctions  as  weighting 
functions  and  integrate  to  get: 

O  *  ^  s^(irr^j')  [c;'(y)  -  iVc;fy)tcfr)KW^.t-3) 

Again  the  only  way  this  equation  can  he-  satisfied  is  if  each 
term  for  constant  i  is  equal  to  zero: 

C,*(y)  t  C;fy)£  (4.4.4) 

Where  F(x,y)  is  Eq  4.1.16.  This  equation  can  be  changed  into 
(see  Eq  4.3.13): 

cay)-  y C;U)+  ]W Gfy'){F(*s,y)-  (/</( 4.4.5) 

•  6fy,/0 

Eq  4.4.4  can  be  turned  into  the  following  matrix  equation: 


q.  y  C;  (i)  +  HOC;  (4.4.6) 

”  F(x0^)GfyK,ye)  (4-4'7) 


Eq  4.4.6  can  be  solved  to  get: 


(l  -  N  D)Cj  -  y  CjO)  (4.4.8) 

This  equation  can  be  solved  by  the  Gauss  elimination  method. 

4 . 5  Results 

This  section  presents  the  results  of  solving  of  Eq 
4.1.15  by  the  Method  of  Central  Finite  Diffferences  and  the 
Method  of  Weighted  Residuals  using  a  polynomial  series  ( Eq 
4.3.1)  and  sine  series  (Eq  4.4.1)  to  reduce  the  partial 
differential  equation  to  an  ordinary  differential  equation 
which  is  then  solved  by  the  Green's  Function  method.  Table  7 
shows  the  results  of  the  CFD  method  and  shows  the  exact 
result  on  the  line  below  for  2x2  and  3x3  interior  nodal 
points.  Table  8  shows  how  the  relative  percentage  error  near 
the  center  decreases  as  the  number  of  nodal  points 
increases.  6x6  nodal  points  was  as  large  as  the  program  run 
on  this  Apple  III  could  handle.  The  CFD  method  generates  a 
tridiagonal  matrix  of  size  N*N  by  N*N  which  needs  a  number 
of  operation  of  order  N*N  to  solve  on  a  computer. 


Table  7.  2-D  CFD/Exact  Solutions  at  Each 
Point  for  N  =  2  and  N  =  3 


x=. 3333 

x= . 6666 

y  = .  6666 

0.169342 

0.166667 

0.196272 

0.  192450 

y=. 3333 

0.483753 

0.481125 

0.561300 

0.555556 

x=. 2500 

x= . 5000 

x=. 7500 

y= . 7500 

0.161330 

0.160041 

0.247200 

0.244981 

0.201088 
0. 198874 

y=. 5000 

0.082402 

0.081660 

0.126317 

0.125000 

0.102830 
0. 101474 

y= . 2500 

0.022097 

0.022097 

0.033879 

0.033825 

0.027586 

0.027459 

Table  8.  Relative  Percentage  Errors  Near 
the  Center  for  2-D  CFD  for 
Various  N 


N  Error(%) 


2  1.98 

3  1.05 

4  0.78 

5  0.50 

6  0.40 


Table  9  shows  the  results  using  the  polynomial  series,  the 
sine  series,  and  the  exact  result  (in  that  order)  at  each 
nodal  point  for  2x2  and  3x3  interior  nodal  points.  Table  10 
shows  the  examples  of  the  best  and  worst  relative  percentage 
errors  for  both  methods  for  2x2,  3x3,  4x4,  and  6x6  interior 
nodal  points . 

As  can  seen  the  MWR  method  combined  with  the  Green's 


Function  method  did  not  give  very  good  results.  In  fact,  the 
best  results  came  in  the  2x2  case  and  got  worse  for  for 
larger  results.  Besides  the  bad  accuracy,  another  problem  is 
that  the  MWR  plus  Green's  Function  method  uses  considerably 
more  operations.  The  differential  equations  for  c ( y )  (Eq 
4.3.6  and  Eq  4.4.3)  depend  on  the  value  of  the  x:  chosen  in 
the  Dirac  Delta  weighting  function  used  in  the  Collocation 


Method  ,  so  there  are  N  equations  for  N  xj values.  Each  of 
these  equations  generates  a  dense  matrix  of  order  N  which 
require  N*N*N  operations  each  to  solve.  As  a  result  the 
total  number  of  operations  required  is  of  order  N  to  the 
fifth  for  MWR /Co 1 1  oca t ion  plus  Green's  Function  method 


versus  N*N  for  the  CFD  method. 


Table  9.  MWR  Results:x  (1-x),  sin(i/Tx),  and  Exact 
at  Each  Point  for  N=2  and  N  =  3 


x=. 3333 

x  = .  666 7 

.6667 

0.057980 

0.055061 

0.055556 

0.200274 

0.192389 

0.192450 

.  3333 

0.057221 

0.041218 

0.048113 

0.018779 
0. 152943 
0.166667 

x=. 2500 

x=. 5000 

x  =  .  7500 

.7500 

0.028060 

0.032901 

0.027459 

0. 104077 
0. 118357 
0.101474 

0.202004 

0.216602 

0.198874 

.  5000 

0.028303 

0.030039 

0.033825 

0.110059 

0.114432 

0.1.2500 

0.230413 

0.233472 

0.244981 

.  2500 

-0.027734 

0.018811 

0.022097 

-0.077424 

0.072418 

0.081660 

0.027208 

0.150224 

0.160041 

Table  10.  Relative  Percentage  Errors  for  2-D  MWR 

for  Different  N  and  Two  Weighting  Functions 


N  x  (1-x)  sin  (  i  7T  x  ) 

Best  Worst  Best  Worst 


2  4.07  18.93 

3  1.57  -225.53 

4  0.66  -164.09 

6  -0.70  101.90 


-0.03  -14.33 
-8.  5  19.82 
-3.61  51.89 
-4.89  189.95 


4.6  Separation  of  Variables 


The  method  of  separation  of  variables  can  be  used  to 
get  the  analytic  solutions  to  partial  differential 
equations.  Another  possible  method  of  solving  a  2-D  partial 


diffential  equation  would  be  to  use  the  method  of  separation 
of  variables  to  turn  the  partial  differential  equation  into 
two  ordinary  differential  equations  which  are  then  solved 
using  the  Green's  Functon  method.  Consider  a  general 

heat  conduction  problem  given  in  figure  5. 


Figure  5.  General  Heat  Conduction  Problem 

-  (4.6.D 

Boundary  Conditions:  7  (&>/)  =  Tfck>4)*0(4.6.2) 

T6*,/)  =  H-t* 0 

A  partial  differential  equation  can  be  only  separated 
if  it  takes  the  following  form  (Ref  4:498-499): 

U,y)  +  (4.6.3) 


Comparing  Eq  4.6.1  and  Eq  4.6.3  we  see  that  only  if 


F(x,y)  takes  the  following  form  will  Eq  4.6.1  be  separable: 


f  +  f)  T  (x^y) 


Letting  T  be  separable  we  get  the  following  results: 


T(kvy)  9  Wy)  (4.6.5) 

X  l  ( &)  4 '  [  X  (*&)  —  O  (4.6.6) 

Boundary  Conditions:  X(0)*XC()-O  (4.6.7) 

Y"(y)+  ,,.0.8) 


Boundary  Conditions:  (.0)  ~  YCl')  —  C>  (4.6.9) 


Eq  4.6.6  and  Eq  4.6.7  could  now  be  solved  using  the  Green's 
Function  method.  The  problem  that  arises  is  that  there  is  a 
pair  of  differentila  equations  (Eq  4.6.6  and  Eq  4.6.8)  to  be 
solved  for  each  choice  of  the  separation  constant  a*a.  As  an 
example,  the  2-D  problem  that  we  have  been  considering  in 
the  previous  sections  (Eq  4.1.12)  becomes: 


(4.6.10, 

Boundary  Conditions:  )f(0)  JfC ()  ~  C>  (4.6.11) 

yty  +  ly?-**  yZc.*!2)Y(y)  =  6 

Y  *  (4.6.12, 

Boundary  Conditions:  Wtf)  »  Y(()  ~  &  (4.6.13) 


For  this  problem  there  is  only  one  valid  choice  for  the 
separation  constant: 


(4.6.10) 


The  difficulty  that  arises  in  the  general  problem,  is 
finding  the  proper  choice(s)  for  the  separation  constant  and 
the  proper  weighting  of  each  valid  solution  such  that  the 
weighted  sum  of  products  of  individual  solutions  matches  the 
boundary  conditions.  Combining  this  difficulty  with  the 
fact  that  potentially  a  large  number  of  ordinary 
differential  equations  will  have  to  be  solved  we  see  that 


this  method  is  not  a  very  viable  or  attractive  method  of 
solving  partial  differential  equations. 


-  *  ’  a  n  C  < 


5.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  purpose  ot  this  thesis  was  to  investigate  a  number 
of  techniques  that  can  be  used  to  solve  ordinary  and  partial 
differential  equations.  The  main  interest  was  on  the  use  of 
Green's  Functions  in  solving  differential  equations.  The 
Method  of  Weighted  Residuals  was  also  used.  The  method  of 
central  finite  differences  was  used  since  it  is  a  commonly 
used  method  and  it  was  desired  to  see  how  the  other  methods 
compared  to  it. 

5.1  One  Dimensional  Problem 


A  one  dimensional  equation  was  used  that  could  be 
solved  using  both  the  Laplacian  Green's  Function  and  the 
Helmholtz  Green's  Function.  The  Laplacian  solution  is  an 
integral  equation.  The  one  dimensional  equation  also  had 
nice  solutions  using  the  Galerkin  and  Collocation  methods 
which  are  subtypes  of  the  Method  of  Weighted  Residuals.  The 
Galerkin  and  Collocation  methods  both  gave  very  accurate 
results  for  only  a  small  number  of  points  used.  The 
Collocation  method  gave  better  results  than  the  Galerkin 
method.  The  central  finite  difference  method  was  easy  to 
accomplish  for  this  problem.  The  results  were  that  CFD 
method  and  the  Laplacian  Green's  Function  method  gave 
exactly  the  same  results  (to  within  accumulated  round-off 
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error).  This  result  was  surprising  considering  the 
completely  different  matrices  in  each  method.  However,  a 
previous  thesis  had  also  gotten  a  similiar  result.  The 
Helmholtz  Green's  Function  method  gave  slightly  better 
results  than  the  CFD  and  Laplacian  methods.  The  Galerkin  and 
Collocation  methods  both  gave  much  better  results  than  the 
CFD  or  either  Green's  Function  method. 

5.2  Two  Dimensional  Problem 


In  the  two  dimensional  case  a  partial  differential 
equation  was  chosen  that  was  separable  and  had  the  form  of 
the  Helmholtz  Eq .  The  CFD  method  gave  reasonably  good 
results.  Using  the  Method  of  Weighted  Residuals  to  reduce 
the  partial  differential  equation  to  an  ordinary 
differential  equation  gave  very  bad  results  for  both 
approximations  that  were  tried.  The  method  of  separaion  of 
variables  was  considered  as  a  way  of  reducing  a  partial 
differential  equation  to  ordinary  differential  equations 
which  would  then  be  solved  using  Green'  Function  method. 

This  method  was  not  actually  solved  because  of  difficulties 
in  choosing  the  separation  constants  and  in  weighting  the 
various  solutions  that  each  sepration  constant  generates.  A 
two  dimensional  Green's  Function  solution  was  looked  at 
during  thesis  preparation  but  the  results  were  not  presented 
in  this  thesis.  The  main  difficulty  was  in  the  summation  of 


the  infinite  semi- periodic  series  expansion  for  the  Green's 
Function. 

5.3  Recommendations 

When  solving  a  one  dimensional  ordinary  differential 
equation,  the  Weighted  Residuals  Methods  should  be 
considered.  They  appear  to  give  very  good  results.  They  have 
the  advantage  that  an  approximate  functional  form  is 
generated  which  can  be  used  to  generate  an  answer  at  any 
point  and  not  just  at  the  nodal  points  used.  One 
disadvantage  of  Weighted  Residual  Methods  is  that  the  sizes 
of  the  relative  errors  vary  considerably  over  the  range  of 
integration.  The  Method  of  central  finite  differences 
involves  fewer  operations  than  Weighted  Residual  Methods  and 
the  Green's  Function  methods.  For  simple  problems,  it  is 
preferable  to  the  Green's  Function  methods.  For  really 
complicated  problems  the  Green's  Functions  methods  may  give 
better  results,  but  this  was  not  actualy  shown  to  be  true. 

Using  a  Weighted  Residual  method  to  reduce  a  two 
dimensional  partial  differential  equation  to  a  one 
dimensional  differential  equation  which  is  then  solved  by 
Green's  Functions  does  not  appear  to  work.  The  central 
finite  difference  method  is  much  preferable  as  a  method  of 
solution  to  the  two  dimensional  problem. 


MUSS 


D 

.  Kr  a 

jci 

ut 

ion  o 

f  B 

an 

ics  . 

Par 

t  e 

Ha  j  d  i 
Me  t  ho 
Struc 
Equat 

Methods  in  Engineering,  4_:509-522  (  1972). 

Hajdin,  N.  and  D.  Krajcinovic.  "Integral  Equation 
Method  for  solution  of  Boundary  Value  Problems  of 
Structural  Mechanics,  Part  II.  Elliptic  Partial 
Differential  Equations",  International  Journal  for 
Numerical  Methods  in  Engineering . 4 : 523-539  (1972). 

Nelson,  G.  Numerical  Methods  and  the  Solution 
of  Boundary  Value  Problems.  Wright  Patterson  Air  Force 
Base:  AFIT  Thesis,  1979. 

Hildebrand,  F.  Methods  of  Applied 

Mathematics .  New  York:  Prentice-Hall,  Inc.,  1952. 
Morse,  P.  and  H.  Feshbach.  Methods  of 

Theorectical  Physics:  Part  I.  New  York:  McGraw-Hill, 
1953. 

Mathews,  J.  and  R.  Walker.  Mathematical 

Methods  of  Physics  2nd  Ed.  Menlo  Park,  California: 

W . A . Ben jami n ,  Inc.,  1970. 

Jackson,  J.  Classical  Electrodynamics  1st  Ed.  New 
York:John  Wiley  &  Sons,  Inc.,  1962. 


Finlayson,  B.  The  Method  of  WeiRhted 
Residuals  and  Variational  Principles. 
Academic  Press,  1972. 


New  Yor k : 


Martin,  H.  and  G.  Carey,  Introduction  to 
Finite  Element  Analysis:  Theory  and 
Application .  New  York:  McGraw-Hill,  1973. 

Houstis,  E.,  R.  Lynch,  and  J.  Rice.  "Evaluation  of 
Numerical  Methods  for  Elliptic  Partial  Differential 
Equations",  Journal  of  Computational  Physics. 
27:323-350  (1978). 


VITA 


William  H.  Amelng  was  born  on  13  January  1953  in 
Monroe,  Michigan,  the  son  of  Howard  W.  Ameling  and  Dorothy 
A.  Ameling  (Toeppe).  He  was  a  National  Merit  Finalist  and 
graduated  from  Ross  High  School  in  Fremont,  Ohio  in  1971.  He 
attended  Case  Western  Reserve  University  in  Cleveland,  Ohio 
from  1971  to  1975  and  received  a  B.S.  degree  in  Astronomy 
with  Honor.  He  then  attended  University  of  Chicago  in  the 
graduate  physics  program  for  three  years.  In  September  1978 
he  left  University  of  Chicago,  joined  the  US  Air  Force,  and 
went  to  Officer  Training  School.  He  was  commissioned  in 
December  1978.  His  first  assignment  was  to  Headquarters  • 
Space  and  Missile  Systems  Organization  (SAMSO)(now  HQ  Space 
Division)  at  Los  Angeles  Air  Force  Station.  He  served  in  the 
Deputy  for  Technology  (SD/YL)  as  a  project  officer  for  the 
Spacecraft  Charging  At  High  Altitudes  (SCATHA)  Progam  for  a 
year  and  a  half.  He  was  then  transferred  to  the  Missile 
Surveillence  Technology  (Measurements)  Program  as  a  project 
officer  for  two  years.  He  was  assigned  to  the  Air  Force 
Institute  of  Technology  in  June  1982  in  the  Engineering 
Physics  Program. 

Permanent  Address:  619  Ewing  St. 


Fremont,  Ohio  43420 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


IEPORT  SECURITY  CLASSIFICATION 

Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


2b.  declassification/downgraoing  schedule 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

AFIT/nFr/ri!/S3D-l 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  ?cr  ’public  release; 
d  iotriV.it  i  >n  .ml  i m  i  ted  . 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OP  PERFORMING  ORGANIZATION 

Air  Force  Institute 
of  Technolopv 

6b.  OFFICE  SYMBOL 
(If  applicable) 

aftt/fmp 

6c.  AOORESS  (City.  Slate  and  ZIP  Code) 

Uriqht-Patterson  AFB,  on  4 5433 

8a.  NAME  OP  FUNDING/SPONSORING 

8b.  OFFICE  SYMBOL 

ORGANIZATION 

(If  applicable) 

Air  Force  Materials  hah 

AF?TJ/TtLB’« 

Be.  ADDRESS  (City.  State  and  ZIP  Code ) 

Uripht-Patterson  AFB,  aw  45433 

7b.  ADDRESS  (City.  State  and  ZIP  Code) 


10.  SOURCE  OF  FUNDING  NOS. 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

NO. 

12.  PERSONAL  AUTHOR(S) 

Arielinp,  Villian  F. 


^  TYPE  OF  REPORT 

Masters  Thesis 


IS.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 

FROM  H.3-03  TO  33-12 


14.  DATE  OF  REPORT  (Yr..  Mo..  Day) 

ln83  Mecomher  5 


TITS 


JVBJiovod  Iot  public  release:  1AW  APR  190-17. 


FIELD 

GROUP 

12 

01 

20 

13 

COSATI  COOES 


Dean  ler  I  .  r  .h  cnJ  Prclc'siona!  DevalopmMt 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necesiMei^kMUidetsti^y armbtafsk  ty/mber) 

■Hpirtfnl  TMf ^r»rnn t ini  Fruatiorn,  ^rcr-n  *  ’•  uricti.  or  s , 
17  in  it  ft  Flerent  .A<nalvsis,  TTfticrhtin<t.  Functions 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

Several  methods  of  numerically'  snlvirf*.  differential  conations  ir  one  and  tie  dimensions 
vere  comnared.  The  main  method  used  ran  a  Green's  Function  intepral  solution.  'T‘l:is  *-’otv,od 
i*a3  compared  with  a  central  finite  difference  (CVp)  solution  and  **itV  tre  b’p.iphtfid 
Residual  methods:  Oalerhin  and  Collocation.  In  a  one  di men sior.al  problem  the  Laplncier 
Creen's  Function  method  j»nve  results  that  were  the  same  as  the  C^D  method..  The  "elrholts 
Green's  Function  method  pavn  slightly  'otter  results  than  th.e  T.aplncinu  Freon's  Function. 
The  Calerkin  and  Collocation  methods  pave  much  1  otter  results  than  the  rrcen's  Function 
For  a  two  dimensional  problem  the  rpp  method,  pave  pood  results.  Tho  Collocation  method 
was  used  to  reduce  the  oartial  differential  equation  to  an  ordinarv  di pf  emntia.1  eouation 
which  was  then  solved  by  the  Creen's  Function  method.  T’-i-?  method  turned  out  to  not  ’-or1-. 
Separation  of  Variables  to  reduce  the  partial  differential  nmirtion  to  an  ordinarv 
differential  eouation  which  is  then  solved  1 v  t' -n  Freer ' «  Function  method  was  also 
considered,  however,  this  approach  was  shown  not  to  1  r  useful  fc*  mo*t  ti*o  dimensior.nl 
problem. _ 


2*i^HSTRI  0UTION/A  V  A I  LA  81  LI  T  Y  OF  ABSTRACT 
UNlLASSiFIEO/UNLIMITEO  &■  SAME  AS  RPT  □  OTIC  USERS  □ 


22a.  NAME  OF  RESPONSIBLE  INOIVIOUAL 

Fap lan,  Bernard 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


22b  TELEPHONE  NUMBER 
( Include  .4nea  Code ) 


22c  OFFICE  SYMBOL 
A  T  TT  t  pTT> 


DD  FORM  1473,  83  APR 


EOITION  OF  1  JAN  73  IS  OBSOLETE. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


